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^ ■ Abstract 

cn ■ We discuss the phenomenology of strange-quark dynamics in the nucleon, based on exper- 

imental and theoretical results for electroweak form factors and for parton densities. In 
particular, we construct a model for the generalized parton distribution that relates the 
[^^ ' asymmetry s{x) — s{x) between the longitudinal momentum distributions of strange quarks 

^D ■ and antiquarks with the form factor Ff (t) , which describes the distribution of strangeness 

in transverse position space. 

X 

1 Introduction 

The distribution of strange quarks and antiquarks is a nontrivial aspect of nucleon structure. Whereas 
the presence of these non-valence degrees of freedom is not surprising, given that gluons can split 
into ss pairs, their relative abundance compared with uu and dd pairs reflects the role of quark 
masses in nonperturbative dynamics. Furthermore, asymmetries in the distribution of s and s are 
not generated by the simple splitting g ^ ss and hence are footprints of more subtle dynamical 
mechanisms. Quantities that have received considerable attention in the recent literature are form 
factors of electroweak currents, which are accessible through parity violation in elastic lepton-nucleon 
scattering, and the difference between the momentum distributions of strange quarks and antiquarks, 
which has in particular shown to be relevant for the determination of the weak mixing angle from 
deep inelastic neutrino-nucleon scattering [1]. In the present work we point out interrelations between 
different measures of strangeness and connect two of them quantitatively in a particular model. 



A number of quantities related to strangeness in the nucleon are matrix elements of local operators. 
In view of our remarks in the previous paragraph, it is important to note the behavior of these 
operators under charge conjugation. In particular, the electromagnetic current is C odd and hence 
sensitive to the difference between contributions from s and s. In contrast, operators like the axial 
vector current, the energy- momentum tensor or the scalar currenl|l| are C even and thus add the 
contributions from s and s. Large values of nucleon matrix elements would be more surprising for 
C odd operators than for C even ones, since for C odd operators they necessitate important effects 
beyond simple g ^> ss fluctuations. 

The unpolarized parton densities s(x), s{x) and their longitudinally polarized counterparts As(x), 
As(x) are expectation values of nonlocal operators and give the momentum distribution of strange 
quarks or antiquarks in a fast moving nucleon. Specific moments of these distributions are associated 
with local operators of definite C parity, as we will specify shortly. 

A suitable framework to discuss relations between various quantities describing nucleon structure 
is provided by generalized parton distributions. They are matrix elements of the same nonlocal 
operators that define the usual parton densities, but taken between proton states of different momenta. 
Throughout this work we consider these distributions at zero skewness ^ = 0, and for brevity we will 
not display this variable. For our discussion it is useful to introduce distributions 

H^x,t) = -Hi{-x,t), E'>{x,t) = -E'i{-x,t), H^{x,t)=Hi{-x,t), (1) 

where the different signs on the r.h.s. reflect the different behavior of vector and axial vector operators 
under charge conjugation. H''{x,t), E'^{x,t) and H''{x,t) respectively correspond to H'^{x,^ = 0, i), 
E'^(x,^ = 0,t) and H''{x,^ = 0, t) in the notation of [3,4]. Taking t = and x > we obtain the 
usual quark and antiquark densities of the proton as 

H'^{x, 0) = q{x) , H^{x, 0) = q{x) , H^{x, 0) = Aq{x) , H^{x, 0) = Aq{x) . (2) 

A two-dimensional Fourier transform with respect to t gives the so-called impact parameter densities 

d^A _.^.^ _. . ... _. .. f d^A 



q{x, b) = J ^ e-"-^ H'^ix, t = - A^) , q{x, b) = J j^ e^^'^ H^{x, t = - A^) , (3) 

which specify the spatial distribution of quarks or antiquarks with longitudinal momentum fraction 
X in the transverse plane, where the impact parameter b is the transverse distance of the parton from 
the center of the proton [5]. Impact parameter densities Aq{x, b), Aq(x, b) for longitudinally polarized 
quarks and antiquarks in a longitudinally polarized proton are obtained from H'^{x, t), H^{x, t) in full 
analogy to ^. The Fourier transforms of E'^{x,t), E^(x,t) describe the dependence of the impact 
parameter distribution of unpolarized quarks or antiquarks on transverse nucleon polarization [5]. 

The distributions just introduced are connected with the form factors mentioned above by sum 
rules, i.e. by integrals over the momentum fraction x. In particular, the lowest moments 

Fl{t) = I dxH'{x,t) = f dx [H%x,t) - H%x,t)] , (4) 

J-i Jo 

F|(t) = / dxE\x,t) = I dx [E'{x,t) - E'{x,t)\ (5) 

give the strange Dirac and Pauli form factors, which are defined as 

{p{p')\srs\p{p)) = Fm^p'h^uip) + Fi{t)n{p') ^^^^^b), (6) 



^We recall that the scalar current ss is relevant in connection with the pion-nucleon a term, see e.g. [2]. 



where t = {p — p') and mp is the proton mass. Then' norniahzation is 

Ff(0) = 0, Fi{0) = fis, (7) 

where the first condition reflects that the proton has no net strangeness, whereas the second condition 
involves the strangeness magnetic moment fig- Note that the contributions of Ff and F^ to the 
electromagnetic form factors of proton and neutron appear with a charge factor e^ = —1/3. The sum 
rules ^ and ^ involve the difference of quark and antiquark distributions, as it must be for form 
factors of the current s^^s. Taking the Fourier transform as in ([3]) we see that 

describes the difference between the transverse spatial distributions of strange quarks and antiquarks, 
averaged over their momentum fraction x. Similarly, the Fourier transform of F^ describes the 
different dependence of the impact parameter distributions on transverse nucleon polarization. 
Further important moments are 

^2,0 W = / dxxH'{x,t)= I dxx[H'{x,t)+H'{x,t)], (9) 

which is a form factor of the energy-momentum tensor for strange quarks, and the strange quark 
contribution to the axial form factor, 

/•I ^ /-i ~ 

F|(i) = / dxH%x,t)= dx[H'{x,t)+H''{x,t)], (10) 

J-i Jo 

which contributes to elastic lepton-nucleon scattering via Z exchange. Both form factors belong to 
charge-conjugation even currents and thus sum over quark and antiquark distributions. 

Using ^ we can connect the values of form factors at t = with moments of the usual quark and 
antiquark densities. In particular, the first condition in ([7]) is equivalent with {s — s) = 0, where we 
introduced the shorthand notation 

(/)= / dxfix). (11) 



In contrast, we have nonzero values for 

^2,o(0) = {x{s + s)), FXm = {As + As) , (12) 

which respectively give the fractional contributions of strange quarks and antiquarks to the longitudi- 
nal momentum and to the spin of the proton. There is no analogous relation for the second condition 
in ([7]) since E^{x,t) and E'^{x,t) do not reduce to any parton density for t = 0. Let us however 
mention that their Fourier transforms with respect to t satisfy positivity constraints involving the 
unpolarized and polarized quark or antiquark densities [6]. 

In order to obtain a quantitative feeling for the role of strange quarks and antiquarks in the 
proton, we briefly review in Sect. [2]the current experimental knowledge of the form factors F^ and F| 
and of strange parton distributions. In Sect. [3] we mention a number of approaches to calculate the 
form factors and the momentum asymmetry s{x) — s{x) theoretically, which will indicate dynamical 
mechanisms that can give rise to these C odd quantities. In Sect. H] we develop a model for H^{x, t) — 
H^{x,t) and use it to connect at a quantitative level the asymmetry s{x) — s{x) with the form 
factor Fl{t). According to ([8]) we thus connect the asymmetry between strange quark and antiquark 
distributions in longitudinal momentum with the one in transverse spatial position. Our results are 
summarized in Sect. El 



Table 1: Data for the strange form factors at low —t. 
added in quadrature. We quote results for C^j or G|; + 



Statistical and systematic errors have been 
rjGlf and the equivalent ones for Ff + r]'F2. 



experiment 


-t[GeV2] 






Ft,F^ 




SAMPLE [7] 


0.100 


G%i = 


0.37(34) 


Ff + Fi = 


0.37(34) 


A4 [8] 


0.23 


G% + 0.225 G%j = 


0.039(34) 


Ff + 0.130 F| = 


0.032(28) 


HAPPEX [9] 


0.477 


G| + 0.392 G^^ = 


0.014(22) 


Ff + 0.184 F| = 


0.010(16) 


A4 [10] 


0.108 


G% + 0.10QG%j = 


0.071(36) 


Ff + 0.068 F| = 


0.064(33) 


HAPPEX [11] 


0.091 


G% = - 


-0.038(43) 


Ff- 0.026 F| = 


-0.038(43) 


HAPPEX [12] 


0.099 


G% + 0.080 GIj = 


0.030(28) 


Ff + 0.048 F| = 


0.028(26) 


HAPPEX [14] 


0.077 


G% = 


0.002(16) 


Ff - 0.022 F| = 


0.002(16) 


HAPPEX [14] 


0.109 


G% + omoGij = 


0.007(13) 


Ff + 0.054 F| = 


0.006(12) 



2 Experimental results for strange form factors and parton densities 

2.1 Electromagnetic form factors 

The strange form factors can be extracted from parity violation in elastic electron scattering on a 
nucleon [7-14]. Experiments typically measure a linear combination of the electric and magnetic form 
factors G|; and G^^ at a low value of the momentum transfer. This can of course be converted into 
a linear combination of Ff and F|, using the relations 



G% + tG 



M 



1 + r 



Fo^ 



A/f — (-f r 



'M 



1 + r 



(13) 



where r = —t/Ami. Recent experimental results at low —t are compiled in Table [H Inspection of the 
table reveals that only the strange Dirac form factor is fairly well determined, whereas F| suffers from 
large uncertainties. It is also evident that the recent HAPPEX data [14] have significantly smaller 
errors than the other measurements. Unfortunately the two form factor combinations given in [14] 
are for different values of t. The determination of the individual form factors therefore requires an 
assumption about their t dependence. A simple way to proceed is to ignore the difference in t of the 
two measurements. Prom the corresponding two entries in Table [U one then obtains 



Fl{t ~ -0.1 GeV^) = 0.003(12) 



F^{t ~ -0.1 GeV^) = 0.05(26) . 



(14) 



This result is graphically represented in Fig. [TJ The use of all data near t = —0.1 GeV in Tabled] 
does practically not change F^, whereas F2 becomes substantially larger but stays within the error 
quoted in (fT^ . 

An alternative method has been used in [15], where a parameterization of the small t behavior 
of the strange form factors was fitted to all data below —t = 0.3 GeV . With the results updated 
in [16], the authors of this study obtain [17] 



Fl{t) = -tx 0.02(11) GeV" 



F|(t) = -0.01(25), 



(15) 
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Figure 1: Results from [14] for the strange form factors at t ~ —0.1 GeV^. The dark (violet) band 
represents the result for G^^, and the light (yellow) band the one for C^ + 0.09 G^j. If one neglects 
the difference of the associated t values one obtains the central value given in ()14p , which is shown as 
a bullet. Also shown are the corresponding la and 2a ellipses. 



for —t < 0.3 GeV , which at — t = 0.1 GeV is fully compatible with our simple estimate (|14p. The 
analysis in [15, 16] includes the data of the GO Collaboration [13] with their very fine binning in t, 
which we have not listed in Table [TJ 

We remark that the experimental values quoted here are subject to theoretical uncertainties due 
to the effects of two-photon and of 7Z exchange, which have been discussed in [18-20] and may not 
be negligible. 



2.2 Unpolarized parton densities 

The determination of parton densities (PDFs) from unpolarized hard scattering processes has made 
significant progress in the recent decade, in particular thanks to data with high precision and a large 
kinematical reach from HERA and from the Tevatron. The knowledge of strange distributions is 
much less advanced, because the observables that dominate global fits of parton densities have little 
sensitivity to s or s. This holds in particular for inclusive deep inelastic scattering (DIS) in kinematics 
where photon exchange is dominant. More specific constraints on s and s distributions come from 
fixed-target DIS experiments with z^ and z^ beams. Thanks to such measurements, there have recently 
been dedicated attempts to determine s{x) and s{x) without strong assumptions on their relation 
with the light sea quark distributions u{x) and d{x). 

In Table [2] we give the moments {x{s + s)) obtained in recent PDF extractions. The study in [21] 
was dedicated to strangeness and explored a number of fits at NLO in a^. The extractions in [22,23] 
were performed at NNLO. The table also gives the moment {x{u + d)), whose values range from 
6.21 X 10^^ to 6.79 X 10^^ and thus show a much smaller spread than for {x{s + s)). The ratio of 
{x{s + s)) and {x{u + d)) is between 0.36 and 0.72 and quantifies the suppression of strangeness in 
the light quark sea. We furthermore show the fiavor asymmetry {x{u — d)), which like {x{s — s)) 
is not generated by simple g ^> qq splitting and hence requires more subtle dynamics in order to 
be nonzero. The ratio of {x{u — d)) and {x{u + d)) varies between —7% and —14%. The parton 









CTEQ6.5S [21] 








set 
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2 3 4 


-0 


-1 


-2 


{x{s + s)) 


3.35 


2.46 


4.44 2.45 4.30 


3.72 


3.48 


4.04 


{x{u + d)) 


6.55 


6.79 


6.21 6.74 6.35 


6.54 


6.65 


6.37 


{x{u — d)) 


-0.72 


-0.75 


-0.69 -0.74 -0.69 


-0.74 


-0.76 


-0.45 



Table 2: Results of different PDF fits for the 
moments {x{s + s)), {x{u + d)) and {x{u — d)). 
All numbers are given in units of 10~^ and refer 
to the scale ^ = 2 GeV in the MS scheme. 





Alekhin 06 [22] MRST 2006 [23] 


{xis + s)) 
{x{u + d)) 
{x{u — d)) 


3.40 3.89 

6.56 6.77 

-0.56 -0.92 



densities corresponding to the entries in the table are plotted in Fig. [2j We note that there are no 
experimental constraints on u{x) — d{x) at small x, so that the large difference between the results 
of [22] and [21,23] for this combination of densities is a consequence of the different functional forms 
assumed in the fits. 

The CTEQ6.5S study [21] also performed fits where s{x) and s{x) were allowed to be different. An 
essential input for constraining this difference are the CCFR and NuTeV data for dimuon production 
in u and u DIS [24]. The parameterization in [21] was chosen such that s{x) — s{x) has precisely one 
zero crossing. The resulting momentum asymmetry at scale // = 2 GeV was found to be 



(x(s-s)) 



2.0 X 10^3 

-0.94 X 10-3 

2.9 X 10-3 



(set -0) 
(set -1) 
(set -2) 



(16) 



where set —0 corresponds to the best fit, whereas sets —1 and —2 were chosen to give the smallest 
and largest values of {x{s — s)), respectively. The ratio of {x{s — s)) and {x{s + s)) in the three 
fits has respective values 5.4%, —2.7% and 7.2%, which is somewhat smaller in size than the ratio of 
{x(u — d)) and {x{u + d)) in the non-strange sector. The left panel in Fig. [3] shows the asymmetry 
s{x) — s{x) obtained in these fits. We note that the best fit (set —0) is quite similar to preliminary 
results obtained by the MSTW collaboration [25]. It should be emphasized that a wider range of 
shapes is obtained if one allows for a variation of the small-x behavior of s{x) — s{x), which is not 
well constrained by data. This is documented in a previous study by CTEQ [26] and shown in the 
right panel of Fig. [3l 



2.3 Polarized parton densities 

The polarization of strange quarks and antiquarks in the proton is not well known at present, for 
similar reasons as in the case of unpolarized parton densities. Many determinations of polarized PDFs 
in the literature, such as those in [27-29], are restricted to the structure functions for inclusive DIS 
with electron or muon beams, which does not permit a separate determination of strange densities. 
This is highlighted in the "valence" scenario of [27], which assumes As(j;) = As(x) = at the starting 
scale /i = 0.63 GeV of evolution. This gives a tiny moment (As + As) ~ —4 x 10-^ at ^ = 1 GeV. 
The study in [30] additionally fits RHIC data for vr" production, which has no particular sensitivity 
to strangeness either. A process that is specifically sensitive to strangeness is semi-inclusive kaon 
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Figure 2: Different combinations of parton densities at /x = 2 GeV. Left: comparison of the different 
sets from CTEQ6.5S [21]. Right: comparison of set from CTEQ6.5S [21] with Alekhin 06 [22] and 
MRST 2006 [23]. 
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Figure 3: The strangeness asymmetry distribution at scale fi = 2 GeV. The left panel shows the fits 
of the CTEQ6.5S analysis [21] and the right panel those of CTEQ6 [26]. Sets B and C correspond to 
different assumptions on the small-x behavior of s(x) — s{x), whereas sets —1 and B— (—2 and B-|-) 
have been chosen to minimize (maximize) the moment {x(s — s)). 

production in DIS, which has been measured by HERMES [31]. The study by de Florian et al. [32] 
includes these data and gives two sets of fits corresponding to different fragmentation functions. It 
does not assume flavor SU(3) symmetry in the polarized sea, allowing As(x) to differ from Au{x) and 
Ad{x). We note that a recent analysis of semi-inclusive hadron production by COMPASS reported 
evidence that the polarized sea is not flavor symmetric and that (An) and (Ad) may have opposite 
sign [34]. All analyses performed so far assume As{x) = As(x). 

In Table [3] we list the values obtained in recent analyses for the first moment (As + As). Note 
that 1/2 times this moment gives the contribution of strange quarks and antiquarks to the total spin 
1/2 of the nucleon. The values from the analyses [27-30] have been obtained with parameterizations 
where the polarized sea quark densities are equal, Au{x) = Ad{x) = As(x), so that they should not 
be regarded as specific determinations of the polarization of the strange sea. Rather, they indicate 
that the contribution of sea quarks to the nucleon spin is negative and of moderate magnitude. 
The different results of the study [32] illustrate that a flavor decomposition of this contribution is 
currently affected with considerable uncertainties. The numbers do not suggest that the strangeness 
contribution to the nucleon spin is very much suppressed compared with the light flavors u and d, but 
further data and analyses are clearly necessary to settle this issue. As a further word of caution we 
remark that an important fraction of the moments in Table [3] comes from the region of small x, where 
the polarized densities are not constrained by data. Quantitative discussions are given in [28, 31] . 



3 Theoretical approaches 

3.1 Electromagnetic form factors 

The strangeness contributions to the electromagnetic and axial form factors of the nucleon have 
been studied in a large number of theoretical approaches (with many studies focusing on the strange 
magnetic moment or the electric charge radius) . Detailed reviews and discussions can be found 
in [7,35,36]. In Table U] we list a small number of recent results for the strange Dirac form factor at 



set 


GRSV 2000 [27] 
"standard" 


BB [28] 
3 4 


LSS 05 [29] 
1 2 


AAC 06 [30] 
1 2 


(As + As) = 


: 2(Ag-) 


-0.126 


-0.148 -0.143 


-0.122 -0.140 


-0.10 -0.12 



Table 3: Lowest moments of polarized parton densities. 
All analyses shown set As(x) = As(a;), and all except 
for [32] take Au{x) = Ad{x) = As(x). The correspond- 
ing PDFs have been determined at NLO accuracy in the 
MS scheme and refer to the scale /x = 1 GeV. 



set 


DNS [32,33] 
KRE KKP 


(As + As) 
(An) 
{Ad) 


-0.095 -0.090 
-0.046 0.076 
-0.048 -0.101 



Table 4: Theoretical results for the strange Dirac form factor at to = —0.1 GeV . The numbers for 
Refs. [43-47] have been read off from graphs. The value in the last row has been obtained from G^^{to) 
in [49] and G^/(0) in [48] using the approximation G^/(to) ~ G^j(O), which was also made in Fig. 2 
of [16]. Taking into account that |G^j(to)| < |G^,j(0)j would increase the value of -Ff(to)- 



Approach 




reference 


Ff{t = -0.1 GeV^) 


Perturbative chiral quark model 




[43] 


0.003 


Chiral quark soliton model (vr) 




[44] 


0.063 


Chiral quark soliton model (K) 




[44] 


0.028 


Vector meson dominance 




[45] 


-0.07 


Vector meson dominance 




[46] 


0.014 


Lattice 




[47] 


0.015(5) 


Lattice + measured magnetic moments and charge 


radii 


[48,49] 


0.000(6) 



—t = 0.1 GeV . We find a substantial spread between these results and remark that several of them 
are outside the range —0.009 < Ff < 0.015 obtained from the experimental values (J14p and (jlSp . 

The calculation of strange form factors is challenging in many theoretical approaches. A large 
number of studies are based on the meson cloud picture, where the nucleon fluctuates into a K and 
a S or A. The coupling to the strangeness current then proceeds through valence degrees of freedom, 
namely the s in the kaon and the s in the hyperon. Concerns have been raised about the quantitative 
reliability of such calculations, based on the possible importance of unitarity corrections [37] and of 
higher-mass states [38,39] such as K* mesons [39]. There seems to be no consensus about these issues 
in the literature, see [40,41] and [36]. We will not quantitatively use the meson cloud picture in the 
present work, but use it as a qualitative guide in Sect. 14.21 Chiral perturbation theory provides a 
systematic framework for calculations in terms of hadron degrees of freedom, but its predictive power 
for strange form factors is limited, as discussed in [42]. 

Kaons also play an essential role in chiral quark models such as the one in [43] , where the nucleon 
is described in terms of three constituent quarks coupling to the pseudo-Goldstone bosons. In this 



approach, nonzero strange form factors are due to the splitting of a u or d quark into a kaon and an s 
quark. The chiral quark sohton model [44] does not rely on the constituent quark picture, containing 
as degrees of freedom both quarks and antiquarks coupling to pions and kaons. 

A different approach is based on dispersion relations, which represent the form factors for spacelike 
t in terms of an integral over their imaginary parts in the timelike region. The assumption that the 
dispersion integral is dominated by single vector meson states leads to the vector meson dominance 
approximation, which underlies many calculations of the strange form factors. A typical procedure is 
to fix the relevant nucleon-meson coupling constants from the isoscalar electromagnetic form factors of 
the nucleon and then to predict the form factors of the strangeness current. Such analyses often obtain 
rather large couplings of the nucleon to the u;(782) and the 0(1020), see for instance [45]. These large 
couplings are in strong conflict with determinations from nucleon-nucleon potential studies [50,51] or 
from dispersion relations for forward nucleon-nucleon scattering [52], with SU(6) symmetry [53], and 
in the case of the (p with the OZI rule. To understand why fits of nucleon form factors with a small 
number of vector meson resonances can give large couplings, we consider the simplified case of just 
two mesons with masses mi and m2, 

ai ^ ^^ ^ aiiml-ml)-{ai+a2){t-ml) 
ml —t ml — t mf 771-2 - ("'-i + w-l) * + *^ 

In order to obtain a 1/t^ behavior of the isoscalar form factor at large —t, one must have ai ~ —02 to 
keep the term with t in the numerator small. At t = the form factor is then approximately given by 
ai{m^ — 771^)/ (771^7712), and the small mass difference between (^(1020) and aj(782) forces the couplings 
ai and 02 to be large. Taking into account higher-mass resonances significantly reduces this trend. 
As an illustration, we have fitted the isoscalar nucleon form factors to a sum of contributions from 
(/>(1020), aj(782) and ix'(1420), with or without an additional contribution from w(1650). We require an 
asymptotic behavior Ff + F^ ^ \/t^ at large — t, which provides a linear relation between the different 
meson couplings in generalization of the simple case we just discussed. With the normalization 
condition -Ff (0) + i^f'(O) = 1 this leaves two free parameters if the a-'(1650) is included and a single 
parameter if this resonance is omitted. For the meson- nucleon couplings relevant to the Dirac form 
factors we obtain 

5^iViV = -9.13, 5;;;^^ = 20.6 (18) 

in the fit without w(1650), where the couplings refer to the ground state mesons and are denoted by 
G^j^. and G]^^^ in [51]. Including the a;(1650), we obtain a good description of the data with 

g^^^ = 4m, g^^^ = n.9, (19) 

where g^j^iN ^^ fixed to a value as small as the data permits, in order to minimize the tension with 
the still lower values obtained in [50,52,53]. Taking both couplings as free fit parameters we find 
9sNN ~ ~0-06 and g^j^jj^ = 14.9. Similar values have been obtained in [54]. 

This simple exercise suggests to take with great care the corresponding predictions for strange 
form factors, where contributions from <j) resonances are strongly enhanced compared with those from 
to states. A more realistic treatment requires the inclusion of continuum states, as has for instance 
been done in [55,56]. The analysis in [56] found that the inclusion of the KK and pir continua makes 
the interpretation of a residual (p resonance contribution ambiguous. 

3.2 The strangeness asymmetry s(x) — s{x) 

The meson cloud picture naturally induces an asymmetry in the momentum distribution of strange 
quarks and antiquarks, which was first observed in [57] and underlies many calculations [58-64]. In 

10 



this picture, the densities of s and s in the nucleon are given as convolutions of the longitudinal 
momentum distributions of the kaon or hyperon within the nucleon and the valence distribution of 
s in the kaon or of s in the hyperon. A similar mechanism is realized in chiral quark models [65], 
where the constituent quarks of the proton can fluctuate into a kaon and a strange quark. There is 
a significant spread among meson cloud model predictions for the shape of s{x) — s{x), including its 
sign and the number of zero crossings, see e.g. the comparative study in [63]. The inclusion of K* 
fluctuations in [64] also had a significant effect, changing even the sign of the momentum asymmetry 
{x{s — s)) compared to the result with kaon fluctuations alone. We note that the predictions for 
s — s in such models typically have a zero at a value of x much larger than 0.01 and are thus rather 
different from the results obtained in the PDF fits [21,25]. 

The study [66] pointed out that in perturbative evolution at three-loop accuracy and beyond, 
graphs with three-gluon exchange in the t-channel generate an s — s asymmetry. Starting with 
s(x) = s{x) at the low scale // = 0.51 GeV, the authors of [66] find that s(x) — s{x) for // > 2 GeV is 
positive at small and negative at intermediate to large x, with (x(s — s)) ~ —5 x 10~^. This is much 
smaller than the central fit results in [21,25], which suggests that this perturbative mechanism plays 
only a minor role in the generation of the momentum asymmetry. 

4 Relating the strange Dirac form factor to s{x) — s{x) 

We now formulate a model ansatz for the C odd part of the generalized parton distribution H^ at 
zero skewness, which will allow us to calculate the Dirac form factor Ff from the phenomenologically 
extracted asymmetry s — s of momentum distributions. Following previous studies of generalized 
parton distributions in the non-strange sector [5,68-70], we assume an exponential t dependence with 
an X dependent slope and set 

H%x, t) - H%x, t) = [s{x) - s{x)] e*-^^(^) , (20) 

where for the slope we take the simple form 

/,(x) = a'(l-x)log-, (21) 

X 

which was already proposed in [5]. With ([3]) it is easy to see that the profile function /s(x) has a 
simple physical interpretation in terms of the average squared impact parameter 

, „ f£bb'^\s(x,b)-s(x,b)] 

J d^b [s{x,b) — s{x,b)\ 

associated with the difference between s and s distributions. As shown in [71], a finite average 
transverse size of parton configurations with x ^ 1 in the nucleon requires {b'^)x to vanish at least 
like (1 — x)^ in this limit, which is obviously satisfied for the ansatz (j2ip . 

In the opposite limit x — > 0, we use simple Regge phenomenology as a guide for our parame- 
terization. The form (j2ip corresponds to the behavior H^ — H^ ~ x""^*-', which arises from the 
exchange of a single Regge pole with a linear trajectory a(t) = a(0) + a't, or from the superposition 
of several Regge poles with the same value of a' . This is a generalization to finite f of a small-x 
behavior x~° for the usual parton densities, which is consistent with phenomenology. The leading 
Regge trajectory that can contribute to H'^ — H^ is the one for the (p mesons, and assuming a linear 
form a^{t) = 0^(0) + a't one obtains 0^(0) = 0.13 and a' = 0.84 GeV~ from the masses and spins 
of 0(1020) and 03(1850). This value of a' is close to the one for other meson trajectories, such as the 
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ones for the p and w. The (f) trajectory contributes to soft hadronic scattering processes hke kaon- 
nucleon scattering or photoproduction of the cj) meson. It is however neglected in most analyses of 
these processes (which is justified if the (/)-nucIeon coupling is sufficiently small, see our discussion in 
Sect. 13. ip . An exception is the analysis of the total kaon-nucleon cross sections performed by Barger 
and Olsson [67], who found an intercept 00(0) = 0.33 =b 0.06 of similar size as the result obtained 
from the hadronic spectrum. We emphasize that in our approach we do not need an explicit value 
for the (/>-nucleon coupling, since the normalization in (I20p is fixed by the difference s — s of parton 
distributions. 

We note that the CTEQ6.5S densities at fj, = 2 GeV are well approximated by 

x{s - s) « ax°-2^ (23) 

in the region 10~^ < x < 10~^, with a = -0.031, +0.023, -0.044 for the respective sets -0,-1,-2. 
This corresponds to a(0) = 0.72, which is quite far from the values we estimate for the (j) trajectory. 
There are however no experimental constraints for the behavior of s — s at very small x, and a value 
of a(0) between 0.1 and 0.4 is within the range for which a good description of all relevant data has 
been obtained in the CTEQ study [26]. 

Using an ansatz as in (j20l) for the valence combinations H^ — H^ and H — H , we obtained 
in [69] a good description of the electromagnetic Dirac form factors of proton and neutron. Given the 
wealth of data in this case, we chose in that study more complicated forms than ()2ip for the profile 
functions /u(x) and fd{x)- We find that for 10^^ < x < 0.1 they are both very well approximated by 
the form (j2ip with a' = 1 GeV~^, which remains close to fu{x) for x > 0.1. Given the fast decrease 
of s{x) — s{x) with X, the small-x region turns out to be most important for our calculation of Ff. 

We take the ansatz ()20p with the CTEQ6.5S densities [21] at fj, = 2 GeV as input, where the chosen 
scale is to be considered as a compromise between a small value appropriate for arguments based on 
hadronic Regge phenomenology and a large value where the densities are sufficiently constrained by 
experimental data. In the left panel of Fig. [4] we show the values of Fl{t) obtained with the best fit 
(set —0) and with the alternative fits (sets —1 and —2). The central curves are for a' = 1 GeV~ in 
(I2ip . and the bands correspond to a' between 0.85 GeV~ and 1.15 GeV^ . We regard this variation 
as an estimate of the parametric uncertainty within our model, with the lower value corresponding to 
the estimate of the (j) trajectory from the meson masses. In the following we refer to the result with 
q' = 1 GeV~ and CTEQ6.5S set —0 as our default prediction. 

If instead of taking (j2ip we set fs{x) equal to the profile functions fu{x) or fd{s) obtained in [69], 
the form factor lies within the bands in the figure, except in the region where |-Ff (t)| has its maximum. 
In that region, the difference of Ff{t) obtained with the different profile functions just mentioned is at 
most 5%. As a further alternative, we have made the ansatz ()20p at scale // = 1.3 GeV, which is the 
starting scale of the CTEQ parameterizations. Taking the profile function (j2ip with a' = 1 GeV~ , we 
again obtain values within the bands of Fig. [H except for deviations of up to 5% around the maximum 
of l-Ff (t)|. Clearly, the largest spread in predictions for Fl{t) within our model is due to the different 
parton densities used as input. To further explore this, we have taken the parameterizations from the 
CTEQ6 study [26] at fj, = 2 GeV, which provides a wider range of shapes as we have seen in Fig. [3l 
The resulting curves for Ff(t) are shown in the right panel of Fig. HI We recall that sets —1 and —2 
in [21] and sets B— and B+ in [26] were chosen to minimize or maximize the moment {x{s — s)). 
They are hence not preferred, although consistent with the data fitted in [21,26]. 

We see that in all cases the form factor Ff{t) is quite small and fully compatible with the estimates 
()14p and (|15p extracted from experiment. We remark that for most of our curves, a linear behavior 
Ff{t) ex t as in ()15p is not a good approximation for —t much above 0.1 GeV . 

It is instructive to compare our results with another small nucleon form factor, namely the Dirac 
form factor F^ (t) of the neutron. Figure [5] shows data together with the default fit from [69] , which 
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Figure 4: The strange Dirac form factor obtained from the sum rule ^ with the model ansatz in 
(I20D and (EB. Central curves are for a' = 1 GeV'^ and bands for 0.85 GeV~2 < a' < 1.15 GeV-^. 



The corresponding parton densities are (from top to bottom): CTEQ6.5S sets —2, —0, —1 in the left 
panel and CTEQ6 sets B-|-, B, C, B— in the right panel. The data point is our estimate (H 
and the dotted curve corresponds to the parameterization of F^ discussed in the text. 
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Figure 5: Data for the neutron Dirac form factor and fit from [69]. The dotted line corresponds to our 
default prediction for Ff, obtained with a' = 1 GeV^ and CTEQ6.5S set —0. The curve is rescaled 
such that its minimum coincides with the one of -F". 
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we already mentioned in connection with the profile functions fuix) and /d(x). The same parame- 
terization multiplied by —0.5 is shown as a dotted curve in Fig. |H We recall at this point that 

Fr(t) = iFtit) - iFr(t) - iFlit) , (24) 

where the labels on the r.h.s. indicate the quark flavor contributions to the Dirac form factor Ff of 
the proton. The fit in [69] neglected the strange contribution in F^ and in F^. We see in Fig. [J] that 
our estimates for Fl are at most half as large in magnitude as -F" for —t = 0.2 GeV^, so that at that 
point —hFf contributes at most 1/6 to the neutron form factor. For higher —t we find that Ff{t) 
decreases faster than Fi{t), which we will explain shortly. Only at small t do we find a stronger 
influence of the strangeness contribution. If our estimate is correct, this is of relevance for the flavor 
analysis of the Dirac radius of the neutron, which in more familiar terms can be expressed through 
the electric radius and a contribution from the magnetic moment, 



{rYi = ^j/iit) 






t=0 



2mS 



{tYe + 0.127 fm^ . (25) 



Concerning the proton form factor, we find that the values of Ff{t) shown in Fig. U] amount to at 
most 3% of Ff{t) at any t. 

4.1 The shape of Fi"(t) 

Let us now discuss the general features of Ff{t) that emerge with our model ansatz, where we have 

Ff{t) = [ dx [s{x) - s{x)] e*-^^(^) , (26) 

Jo 

j^Flit) = I dx [s{x) - s{x)] fs{x) e*^^(^) . (27) 

For the neutron form factor the situation is slightly more complicated even if we neglect the strangeness 
contribution, since the fit in [69] required different profile functions for u and d quarks. Since their 
difference is only moderate, the discussion of F^{t) is however quite similar. 

With fs{x) being a decreasing function, the factor e^'"^^' increasingly suppresses small x values in 
the integral (I26p when —t becomes larger. For increasing —t the form factor Fl is therefore connected 
with s — s at increasing values of x. With the profile function (j2ip we obtain the Drell-Yan relation 
p = i(l + [5) between the powers describing the asymptotic power laws Ff ~ (— t)~P for t — > — cxd and 
s — s ^ (1 — xf for X — > 1 [69, 71]. The difference s — s oi sea quark distributions falls off faster with 
X than the valence distributions u — u and d — d, which are relevant for F", so that one generally 
expects |-F{*(t)| to decrease faster than \F^{t)\ with —t. As Figs. H] and [5] show, this is indeed the case 
in our model. 

In Fig. [6] we show the integrands of Ff and dFf/dt in (j26p and (j27p . The integrands are multiplied 
with X in the plots, so that with the logarithmic scale for x we obtain the form factor or its derivative 
as the area under the corresponding curve. For t = the integrand of Ff is s{x) — s{x), which gives 
a zero integral because of quantum number constraints. The integrand for the derivative dFf/dt has 
an extra factor fs{x), which enhances small x values relative to larger ones. At i = one therefore 
has dFf/dt < at t = if s — s is negative at small x and positive at large x. This is the case for the 
CTEQ fits [21,26] except for sets —1 and B— . As —t increases, the factor e*-*^'^^ suppresses small x 
values in ([27p . and for sufficiently large —t the derivative dFf /dt has the opposite sign than at t = 0. 
For some t one hence obtains a maximum or minimum of Ff{t). The value of —t where this happens 
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Figure 6: Upper plots: the scaled integrands of Ffit) and dFf{t)/dt in ([26]) and (j27|) . obtained with 
our default prediction. Lower plots: the corresponding scaled integrands for F"(t) and dF]^{t)/dt, 
obtained with the default fit in [69] . The values of t in the right panels are as for the corresponding 
curves on the left. 
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is larger for parameterizations of s — s which have the zero crossing at larger x, as one can check by 
comparing Figs. [3] andSl 

Figure [6] also shows the integrands for F" and dF^ /dt for the default fit in [69]. The discussion 
for the sign of the derivative dF^/dt and the presence of a minimum of -F"(t) proceeds in analogy to 
the case of the strangeness form factor. Since the zero crossing oi ^{d — d) — ^{u — u) occurs at much 
larger x than the one of s — s in the CTEQ6.5S parameterizations, whereas the respective profile 
functions are similar, |F"(t)| assumes its maximum at significantly larger —t than |F^(i)|. 

Let us at this point mention the PDFs extracted in [72]. In contrast to the analyses by CTEQ 
[21,26] and MSTW [25], the strangeness asymmetry x{s — s) in [72] has a zero at 0.4 < x < 0.5 and a 
maximum at x ~ 0.7. Taking this distribution with the same profile functions fsix) explored above, 
we obtain a form factor -Ff (t) with a very flat maximum F^ ~ 0.0025 for 1 GeV^ < — t < 2 GeV'^ and 
a slow decrease with —t. In this case, the bulk of the form factor integral (I26p comes from very large 
X, where we deem our model for the profile function fs{x) associated with sea quarks very insecure. 
Since the study [72] used inclusive cross sections for v and P DIS but no dimuon production data 
to constrain s — s, we do not regard such a scenario as strongly motivated. This example illustrates 
however that within our model framework, strong changes in s — s result in qualitatively different 
forms of Ff{t), which may eventually be ruled out by data. 

The relations between the x dependence of s — s and the t dependence of Ff{t) discussed in this 
subsection follow from the general features of our ansatz in (|20p and (j2ip and will also hold for more 
complicated forms. The neutron form factor F^ (t) and the combination ^ (d — d) — g (n — n) of valence 
distributions, which are both much better known than their strangeness counterparts, provide an 
example where these relations are indeed seen and corroborate our prediction for the behavior of 
Ff{t) with a given form of s{x) — s{x). 

4.2 A modified ansatz 

Our ansatz in (j20p is special in that it assumes a t dependence in the form of a global factor multiplying 
s{x) — s{x). It implies that the difference s{x,b) — s{x, b) of impact parameter densities has a Gaussian 
shape in b and in particular does not change sign for given x. One may wonder whether this ansatz 
is too restrictive. The physical picture of meson cloud models for instance suggests that the typical 
transverse position of s is smaller than for s, since the s originates from a kaon, which due to its 
smaller mass tends to be at larger distances than the hyperon containing the s. If this effect is strong 
enough, one may have a node of s{x, b) — s{x, b) in b. 

It is however important to realize that at /i = 2 GeV, where we formulate our model, the individual 
distributions of s and s are not valence-like as they would be in a model valid at low resolution scale. 
For sets —0, —1 and —2 of the CTEQ6.5S parameterization at /i = 2 GeV we find 

x(s + s) «0.2 2;-°-2 (28) 

in the region 10~^ < x < 10~^, which is to be compared with (j23p . For x < 10^^ the ratio (s— s)/(s+s) 
does not exceed 1% in absolute size. It is natural to assume that the bulk of s and s in that region 
is generated through gluon splitting (7 — > ss as described by perturbative evolution. This mechanism 
does not introduce an asymmetry in the transverse distribution of s and s. When introducing a more 
general ansatz for H'^ — H^ than (120p we should hence make sure that the strong cancellation between 
s and s at small x takes place not only in the forward limit but also at nonzero t. With this in mind, 
we explore a variant of (j20p . given by 

H'{x, t) - H'{x, t) = s{x) e*^=(^') - s{x) e*-^'»(^) (29) 
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Figure 7: Left: the difference in the impact parameter distributions for strange quarks and antiquarks 
at X = 0.2. The distributions are multiphed with 27r6, so that the area under each curve gives 
s(x) — s(x). Right: the integrand (j26p of the strange Dirac form factor at —t = 0.5 GeV . The curves 
are for the ansatz in (j29p and ()30p with ol = \ GeV~ and different values of A and A (given in units 
of GeV~^). For the parton densities we take CTEQ6.5S set —0 at // = 2 GeV. 



with 



/s(x) = a'(l — x) log — h Ax{\ — x)^ , 

X 

fJx) = a' (I - x) log - + Ax(l - xf , 

X 



(30) 



where the prefactor x in front of A and A guarantees the cancellation just discussed, as long as —tA 
and —tA are not too large. In the following we take values 2 GeV~ and 4 GeV~ for either A or A. 
With a' = 1 GeV~ this respectively corresponds to a change of fs{x) or fs{x) by a factor 1.2 and 
1.4 at x = 0.2, which one may view as a typical momentum fraction for s and s in a model at low 
scale, where nonperturbative effects could generate an asymmetric distribution in impact parameter. 
In the left panel of Fig. [7] we show the corresponding impact parameter densities s{x, b) — s{x, b) and 
see that they indeed have nodes in b when A or ^ is equal to 4 GeV~ . 

The form factors obtained with this ansatz are shown in Fig. [8l For ^4 > 0, where s is concentrated 
at larger impact parameters than s, we find that Ff is increased in size but not much changed in shape 
compared with our default prediction with A = A = 0. In contrast, we find that for sufficiently large 
A > the form factor Ff changes sign at some finite —t. We can understand this at the level of the 
form factor integrand: for fs{x) > fs{x) the exponential factors in (129^ give a stronger suppression in 
the first term, so that at large enough x and —t one can have H'^{x, t) < H^{x, t) despite s{x) > s{x). 
This is illustrated in the right panel of Fig. [71 

As discussed above, the meson cloud picture suggests that s has smaller rather than larger typical 
impact parameters than s, so that we do not see a particular physics motivation for our examples with 
^ > 0. However, they show that certain nontrivial correlations between the x and b dependence of the 
s and s distributions can have quite drastic effects on Fl{t), which may be observable in experiments 
with sufficient sensitivity and kinematic coverage. 
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4.3 Other form factors 

With our model (|20p for H^ — H^ we can also evaluate the Mellin moment A^ Q(t) = J_i dxx^H^ix^t), 
which is a form factor of an operator with two covariant derivatives between the strange quark field and 
its conjugate. The factor x^ strongly suppresses small x values, and with the CTEQ parameterizations 
for s — s the integral is dominated by an x region where the integrand has a definite sign. As a 
consequence 1^1 Q(t)| decreases monotonically with —t. Its value at t = is tiny, ranging from — 10~^ 
to 4 X 10~^ for the CTEQ6 and CTEQ6.5S parameterizations at /u = 2 GeV. Exceptions are the 
values -6 x 10"^ and lO'^ for CTEQ6B- and CTEQ6B+, respectively. 

We do not attempt here to model the strangeness contributions to the energy-momentum and 
axial form factors, ^2o(^) ^^^ -^a(^)' "which according to ([9]) and (jlOp are respectively related to 
s + s and As + As. These distributions mix with gluons under evolution, which invalidates simple 
ansatze based on Regge trajectories for their small-x behavior. This also holds at finite t [73]. Since 
even the t dependence of H^ + H^ and H + H is barely constrained at present, we see no clear 
guidance for how to model profile functions of H" + H^ and H'^ + H^ . We expect however that 
these distributions have no zeroes in x, which certainly holds for their values at t = according 
to current PDF parameterizations. We therefore predict the corresponding form factors to decrease 
monotonically in absolute size, with values at t = given by the moments in Tables [2] and [31 



5 Summary 

We have discussed several measures of strangeness in the nucleon. Strange quarks and antiquarks are 
not particularly rare in the proton: their contribution {x{s + s)) to the nucleon momentum is only 
suppressed by about a half compared with the one from light flavor antiquarks, {x{u + d)). Their 
contribution (As + As) to the spin of the proton is not well determined at present, but there are no 
indications that it is very much suppressed compared with (An + Ad). 

More subtle quantities are asymmetries between strange quarks and antiquarks, most notably 
the asymmetry between the parton densities s{x) and s(x), and the strange Dirac and Pauli form 
factors Ff[t) and F^{t). A two dimensional Fourier transform of Ff{t) yields the difference of spatial 
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distributions for s and s in the transverse plane, whereas s{x) — s{x) gives the difference of their 
distribution in longitudinal momentum. The two asymmetries are connected via generalized parton 
distributions at zero skewness, for which we have made a model ansatz in order to explore this 
connection quantitatively. Using as an input different sets of s — s distributions extracted by the 
CTEQ Collaboration, we find values of Fl{t ~ -0.1 GeV^) between -0.006 and 0.012, in good 
agreement with current experimental constraints. Many theoretical analyses of the electromagnetic 
nucleon form factors neglect the strangeness contributions. With our estimates this is at most a 3% 
effect for Ff{t). However, l-Ff (t)| might amount to as much as 1/6 of -Ff'(t) at —t = 0.2 GeV^. For 
higher —t the relative contribution quickly decreases, whereas for lower —t it may even be larger. 

The general features of our model ansatz for generalized parton distributions lead to correlations 
between the x dependence of s — s and the shape of Ff{t). The best fits in the PDF extractions 
[21,25,26] yield forms where s — s is negative for small x and positive for large x. With our ansatz, 
this gives a negative derivative dFf{t)/dt at t = and a maximum of the form factor at some value 
of —t. With a zero crossing of s — s at x between 10~^ and 10~^, we find this maximum at —t 
between 0.2 GeV and 0.4 GeV . Analogous correlations are seen to hold between the combination 
■^{d — d) — -^{u — u) of valence quark distributions and the neutron form factor -Ff'(t), which we take as 
support for our predictions. Finally, a rapid decrease of s — s with x reflects itself in a faster decrease 
of Ff(t) compared with F^{t) for large —t. It will be interesting to confront these predictions with 
future data from parity violating electron-nucleon scattering. 
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